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ABSTRACT 

We present a detailed spectroscopic analysis of the luminous blue variable AG Carinae during the 
last two visual minimum phases of its S-Dor cycle (1985-1990 and 2000-2003). The analysis reveals 
an overabundance of He, N, and Na, and a depletion of H, C, and O, on the surface of AG Car, 
indicating the presence of CNO-processed material. Furthermore, the ratio N/0 is higher on the 
stellar surface than in the nebula. We found that the minimum phases of AG Car are not equal 
to each other, since we derived a noticeable difference between the maximum effective temperature 
achieved during 1985-1990 (22,800 K) and 2000-2001 (17,000 K). Significant differences between 
the wind parameters in these two epochs were also noticed. While the wind terminal velocity was 
300 kms~^ in 1985-1990, it was as low as 105 kms~^ in 2001. The mass-loss rate, however, was lower 
from 1985-1990 (1.5 x 10"^ Mgyr-^) than from 2000-2001 (3.7 x 10"^ Mgyr-^). We found that the 
wind of AG Car is significantly clumped (/ ~ 0.10 — 0.25) and that clumps must be formed deep in 
the wind. We derived a bolometric luminosity of 1.5 x 10^ Lq during both minimum phases which, 
contrary to the common assumption, decreases to 1.0 x 10^ Lq as the star moves towards maximum 
flux in the V band. Assuming that the decrease in the bolometric luminosity of AG Car is due to the 
energy used to expand the outer layers of the star (Lamers 1995), we found that the expanding layers 
contain roughly 0.6 — 2 Mq. Such an amount of mass is an order of magnitude lower than the nebular 
mass around AG Car, but is comparable to the nebular mass found around lower-luminosity LBVs 
and to that of the Little Homunculus of Eta Car. If such a large amount of mass is indeed involved in 
the S Dor-type variability, we speculate that such instability could be a failed Giant Eruption, with 
several solar masses never becoming unbound from the star. 

Subject headings: stars: atmospheres — stars: mass loss — stars: variables: other — supergiants — 
stars: individual (AG Carinae) — stars: rotation 



1. INTRODUCTION 

Luminous Blue Variables (LBVs) are characterized as 
high-luminosity, blue (i.e. hot) massive stars showing 
photometric and spectroscopic variability (Conti 1984; 
Humphreys & Davidson 1994; van Genderen 2001; Clark 
et al. 2005). AG Carinae (=HD 94910, Q!2ooo=10:56:11.6, 
'52ooo=-60:27:12.8), one of the brightest, most famous, 
and most variable LBVs, is considered one of the proto- 
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types of this class, presenting all of the above character- 
istics. 

Intense observing campaigns in the past decades have 
determined many of the characteristics of AG Car. The 
star has photometric, spectroscopic, and polarimetric 
variability on a range of timescales spanning days to 
decades. The optical photometric variability of AG Car 
was analyzed in detail by van Genderen et al. (1988), 
van Genderen et al. (1990), Sterken et al. (1996), and 
van Genderen et al. (1997a). These authors determined 
that AG Car has strong photometric changes of up to 
2.5 mag in the V band on a timescale of about 5-10 
years, which are attributed to the S Dor-type variabil- 
ity (van Genderen 1982; van Genderen et al. 1997b; van 
Genderen 2001), and should not be confused with gi- 
ant eruptions similar to that experienced by Eta Cari- 
nae in the 1840's. Superimposed on these variations, 
0.1 — 0.5 mag changes are present on a timescale of 371 
days (van Genderen et al. 1997a), and photometric mi- 
crovariability (0.01 — 0.02 mag) with a timescale of 10-14 
days has been recorded during visual minimum (van Gen- 
deren et al. 1988, 1990). Whitelock et al. (1983) reported 
the near-infrared JHKL lightcurve of AG Car between 
1974-1982, showing that the amplitude of the light vari- 
ations in the near-infrared is comparable to what is seen 
in the V band. 

AG Car is also a spectroscopic variable (Caputo & 
Viotti 1970; Humphreys 1970; Viotti 1971; Wolf & Stahl 
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1982), and the variability on timescale of years has been 
analyzed in different spectral regions: ultraviolet (Lei- 
therer et al. 1994; Shore et al. 1996), optical (Leitherer 
et al. 1994; Stahl et al. 2001), and near- infrared (Groh 
et al. 2007). During the epochs of minimum in the 
visual lightcurve (hereafter referred to as "minimum"; 
my — 8.1; van Genderen 2001), the star is relatively 
hot and has a WNll spectral type (Smith et al. 1994), 
showing strong He I, H I, and N II lines in emission, 
weak He II 4686 A emission, and Si IV 4088-4116 A 
absorption (Stahl 1986; Smith et al. 1994; Walborn & 
Fitzpatrick 2000; Groh et al. 2006). During the maxi- 
mum epochs of the lightcurve (my 6.0; van Genderen 
2001) the star is cooler, and the spectrum is reminiscent 
of that of extreme A-type hypergiants, with strong emis- 
sion of H I, Fe II, and Ti II lines (Wolf & Stahl 1982; 
Stahl et al. 2001). The transition between both phases 
is characterized by the appearance of peculiar features 
in the spectrum, such as absorption-line splitting, strong 
electron-scattering wings in He I and Fe II lines, and ap- 
parent inverse P-Cygni profiles in He I lines (Leitherer 
et al. 1994; Stahl et al. 2001). The variable and high 
polarization degree measured in AG Car motivated spec- 
tropolarimetric monitoring (Schulte-Ladbeck et al. 1994; 
Leitherer et al. 1994; Schulte-Ladbeck et al. 1997; Davies 
et al. 2005, 2006, 2008). 

The presence of a massive, bipolar circumstellar neb- 
ula around AG Car is a testament that a recent {t < 10^ 
years) phase of high mass loss has occurred. The mor- 
phology and kinematics of the nebula were analyzed by 
Nota et al. (1992) and Nota et al. (1995), who deter- 
mined a dynamical age of 8.5 x 10^ years and a high mass 
of ionized nebular material (~ 4.2 Mq), which is likely 
composed of ejecta from the central star (Lamers et al. 
2001). The nebular abimdances were obtained by Mitra 
& Dufour (1990), de Freitas Pacheco et al. (1992), and 
Smith et al. (1997), showing evidence of moderate nitro- 
gen enrichment. Properties of the circumstellar nebula 
were also studied in the mid- and far-infrared, revealing 
an incredibly high dust mass of ~ 0.25 Mq (Voors et al. 
2000), dust temperature between 76 99 K (Voors et al. 
2000), and the presence of large grains of ~ 1/im in order 
to explain the far-IR excess (Voors et al. 2000; Hyland & 
Robinson 1991). Assuming a normal gas-to-dust ratio of 
100, the total nebular mass of AG Car could be as high 
as ~ 30 Mq (Voors ct al. 2000), which would be of the or- 
der of, or even higher than, the mass of the Homunculus 
nebula aroimd Eta Car obtained by Smith ct al. (2003). 
At the time when the AG Car nebula was ejected, the in- 
terstellar bubble around the central star likely contained 
a negligible amount of material compared to the total 
mass of the nebula (Lamers et al. 2001), implying that 
most of the nebular mass was ejected by the central star. 

While insights into the nature of AG Car have been ob- 
tained through the quantitative works of Leitherer et al. 
(1994) and Stahl et al. (2001), the stellar and wind pa- 
rameters of AG Car at each epoch along the S Dor-type 
variability cycle (hereafter S Dor cycle) are not com- 
pletely constrained. Furthermore, it is currently feasible 
to include important physical proc:esses in the radiative 
transfer codes which have had a key impact on the deter- 
mination of stellar parameters of 0-type (e.g. Crowther 
et al. 2002; Hillier et al. 2003; Bouret et al. 2003, 2005; 



Martins et al. 2005; Puis et al. 2006) and Wolf-Rayet 
stars (e.g. Hillier & Miller 1999; Hamann et al. 2006; 
Crowther 2007; Grafener & Hamann 2008). The inclu- 
sion of full line blanketing due to 100 000s of lines, wind 
clumping and charge-exchange reaction, and the simul- 
taneous treatment of the wind and photosphere in the 
model will considerably affect the emerging model spec- 
trum. Together with the possibility of doing a full syn- 
thesis of the spectrum and thus using more diagnostic 
lines, our analysis should provide a more accurate deter- 
mination of the stellar and wind properties of AG Car. 

This work is the first in a series of papers resulting from 
a long-term multi-wavelength spectroscopic and photo- 
metric monitoring campaign to follow AG Car during the 
last two decades. In a previous paper (Groh et al. 2006), 
we presented quantitative evidence that AG Car had a 
high rotational velocity during the minimum phase of 
2000-2003 and, for the first time, we showed that LBVs 
can be fast rotators. In this work, we present, interpret, 
and discuss the results obtained from observations in the 
last two well-documented minimum phases of AG Car; 
namely, during the years of 1985-1990 and 2000-2003. 
Although it is subjective to define the beginning and end 
of minimum phases for stars which have irregular varia- 
tions, such as AG Car, we analyze in this paper epochs 
when the star was fainter than V = 7. This criterion 
corresponds roughly to epochs when He I lines were rel- 
atively strong and displayed a classical P-Cygni profile, 
and time-dependent effects were minimized. 

This paper is organized as follows. In Sect. 2 we 
present the photometric and spectroscopic observations 
obtained during the last two minimum phases of AG Car, 
while in Sect. 3 we describe the spectroscopic evolu- 
tion during those minima. In Sect. 4 we introduce 
the main characteristics of CMFGEN (Hillier & Miller 
1998), the radiative transfer code used to derive the fun- 
damental parameters of AG Car. Sect. 5 briefly dis- 
cusses the distance of AG Car, and how the parame- 
ters are affected by that quantity. The surface chemical 
abundances, mass-loss rate, wind terminal velocity, stel- 
lar temperature, stellar radius, and stellar luminosity for 
each minimum are presented in Sect. 6, and the inter- 
pretation and discussion can be foimd in Sect. 7. The 
conclusions of this paper, which suggest key changes in 
the current paradigm for interpreting the LBV phase and 
their S-Dor cycles, in particular the minimum phases, are 
summarized in Sect. 8. 

The present work will be followed by two accompany- 
ing papers. In Paper II we will discuss the evolution- 
ary status, current mass, proximity to the Eddington 
limit, and whether the bi-stability mechanism is present 
in AG Car. The evolution of the fundamental parameters 
of AG Car during the full S Dor cycle will be analyzed, 
in detail, in Paper HI. 

2. OBSERVATIONS 

2.1. Photometry 

The optical photometry of AG Car was obtained from 
published data from the Long-term Photometry of Vari- 
ables project (LTPV; Sterken 1983; Manfroid et al. 1991; 
Sterken et al. 1993; Spoon et al. 1994; Manfroid et al. 
1995; Sterken et al. 1995), accessed through the online 



On the nature of AG Car I. Minimum phases 



3 



catalog of Vizier/CDS^. The observations were acquired 
using uvby Stromgren filters on the ESO 50 cm La Silla 
telescope, and the typical individual errors are between 
0.003-0.006 mag. The reader is referred to the aforemen- 
tioned papers for further details about the observational 
setup, data reduction, and calibration. Since the LTPV 
project was not active after 1994, CCD V magnitudes 
from 2000-2003 were obtained from the ASAS-3 project 
(Pojmanski 2002). The ASAS-3 data were averaged in 
bins of 15 days, and only high quality (grade A) data 
were used. Although the ASAS-3 magnitudes were mea- 
sured using large apertures (75"), AG Car is by far the 
dominant source of light in the field. 

For comparative purposes, estimates of the visual mag- 
nitude of AG Car from 1981-2004 were compiled from the 
AAVSO archive (A. A. Heiden 2007, private communica- 
tion). When both AAVSO and ASAS-3 data were avail- 
able at a given epoch, their 15-day averaged magnitudes 
were compatible within ~ 0.1 mag. However, the 15- 
day averaged AAVSO visual lightcurve does not exactly 
match the more precise photoelectric V^-band lightcurve 
presented by van Genderen (2001), in particular during 
maximum. This is not surprising, since the AAVSO sam- 
ple is composed by visual estimates (i.e by eye) from 
different observers, and the uncertainty could be large 
during maximum. We estimate that the 15-day averaged 
AAVSO magnitudes have an uncertainty of 0.1-0.2 mag 
during minimum, and perhaps even more during maxi- 
mum. Note that the AAVSO magnitudes are presented 
in this paper only for the purpose of illustrating the gen- 
eral behavior of AG Car during 1981-2004, and due to 
their relatively low accuracy they were not used to derive 
a F-band flux. 

The near-infrared photometry was obtained using the 
75cm reflector telescope from the South African Astro- 
nomical Observatory (SAAO) and JHKL filters. The 
SAAO photometric system is described by Carter (1990), 
and the instrument used was an Mkll photometer with 
a 36"aperture. The typical error for each individual mea- 
surement is 0.02 mag in the J, H, and K bands, and 0.05 
mag in the L band. In principle, the large aperture used 
could include some contribution from the bipolar nebula 
around AG Car. However, magnitudes obtained with 
diaphragm apertures of 54", 36", 24", and 12" do not 
present significant differences, indicating that the nebu- 
lar contribution is negligible (Whitelock et al. 1983). 

Figure 1 displays the AG Car lightcurve encompassing 
its last two minimum phases measured in the K-band, the 
photoelectric Stromgren y-band photometry obtained by 
the LTPV project, the ASAS-3 V-band observations, and 
the photoelectric lightcurve presented by van Genderen 
(2001). For comparison, 15-day averaged visual esti- 
mates made by the AAVSO observers are also shown. 
Table 1 summarizes the photometry of AG Car used in 
this paper. The epochs from the first column correspond 
to those when spectroscopic data, during minimum, were 
available. The magnitudes were obtained through linear 
interpolation between the closest dates which had avail- 
able photometry. 

Fortunately, the time samplings of both the near-IR 
lightcurve and the optical photometry are suitable to fol- 
low the variability of AG Car on timescales of the order 

® http:/ /vizier. u-strasbg.fr/viz-bin/VizieR 
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Fig. 1. — AG Car photoelectric lightcurve from 1982 to 2004 ob- 
tained in the K-band (filled circles connected by a dashed line), in 
the Stromgren y-band by the LTPV project (crosses; Sterken 1983; 
Manfroid et al. 1991; Sterken et al. 1993; Spoon et al. 1994; Man- 
froid et al. 1995; Sterken et al. 1995), and in the V-band by the 
ASAS-3 project (open circles; Pojmanski 2002). For illustrative 
purposes, the schematic V-band photoelectric lightcurve presented 
by van Genderen (2001) (red solid line) and the 15-day averaged 
visual magnitudes measured by the AAVSO observers (dotted line; 
A. A. Heiden 2007, private communication) are also shown. The 
vertical lines mark epochs when spectroscopic data, during mini- 
mum, were available (see Table 2). 

of months, and match the dates when the spectra were 
taken relatively well (Table 1). Variability on timescales 
shorter than a couple of months is definitely present in 
AG Car (van Genderen 2001), however with an ampli- 
tude of < 0.2 mag (van Genderen et al. 1988), and thus 
do not affect the results obtained in this paper. 

2.2. Spectroscopy 
2.2.1. Mmimum of 1985-1990 

Flux-calibrated ultraviolet spectra of AG Car were ob- 
tained by the International Ultraviolet Explorer (lUE) 
satellite in the region 1150-3120 A and gathered through 
the Space Telescope Science Institute (STSCI) public- 
data archive. The lUE absolute-flux calibration is pre- 
cise to about 5% (degrading in the edges of the spectral 
range), which is adequate for our purposes. AG Car was 
observed in various epochs during the 1985-1990 mini- 
mum with the SWP and LWP high-resolution cameras, 
with R ~ 15,000 in the wavelength range of 1150-3200 
A (Table 2, see also Leitherer et al. 1994; Shore et al. 
1996). 

Optical spectra of AG Car were obtained using differ- 
ent telescopes in the southern hemisphere. The Echelle 
spectrograph mounted on the 1.4m CAT/ESO (Chile) 
telescope was used to collect high-resolution spectra from 
1986 June to 1991 January around Ha, He I 5876 A, and 
He II 4686 A, with R=60,000. Most of these data have 
been described elsewhere (Bandiera et al. 1989; Viotti 
et al. 1991; Leitherer et al. 1994). Lower resolution 
data between 3800-4900 A were obtained using the Im 
CTIO/Yale telescope and were discussed by Walborn & 
Fitzpatrick (2000) . Additional lower-resolution observa- 
tions using the longslit Cassegrain spectrograph on the 
1.52m ESO/La Silla telescope were made on 1990 June 
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TABLE 1 

Photometry of AG Car interpolated to epochs when spectra were available during minimum* 



Epoch 


u 


V 


b 


y 


visual 


V 


J 


H 


K 


L 


ALTPV 


AV ^ 


ANIR ^ 


(0.01) 


(0.01) 


(0.01) 


(0.01) 


(0.1) 


(0.05) 


(0.02) 


(0.02) 


(0.02) 


(0.05) 


(days) 


(days) 


(days) 


looc T..1 -in 

l\)oo Jul 19 


8.94 


8.72 


8.37 


7.96 


8.2 




6.45 


6.13 


5.78 


5.34 


112 




-3 


lyob Jun i / 


9.04 


8.78 


8.41 


7.98 


8.1 




6.47 


6.12 


5.80 


5.35 


— 7 




1 


lyoD Jun io 


9.05 


8.79 


8.41 


7.98 


8.1 




6.47 


6.12 


5.80 


5.35 


—6 




2 


1950 Jun 2o 


9.07 


8.81 


8.43 


8.00 


8.1 




6.47 


6.13 


5.81 


5.36 


— 1 




7 


1987 Jan Uo 


9.03 


8.81 


8.43 


8.03 


8.1 




6.52 


6.19 


5.86 


5.36 


192 




—4 


lyo ( Jun iu 


9.00 


8.81 


8.43 


8.03 


8.2 




6.57 


6.24 


5.93 


5.63 


— 183 




12 


lyoY Jul 24 


8.99 


8.80 


8.43 


8.03 


8.1 




6.57 


6.23 


5.91 


5.50 


— 139 




14 


19o9 Mar 2d 


8.98 


8.83 


8.45 


8.07 


8.2 




6.54 


6.19 


5.87 


5.22 


1 




2 


19o9 Uec 26 


8.79 


8.63 


8.26 


7.85 


8.1 




6.38 


6.05 


5.71 


5.25 


1 




—8 


1990 Apr SO 


8.84 


8.72 


8.36 


7.94 


7.9 




6.42 


6.11 


5.81 


5.37 


—64 




—15 


1990 Jun lb 


8.88 


8.72 


8.35 


7.93 


8.0 




6.41 


6.08 


5.77 


5.32 


— 17 




8 


1990 Jun lo 


8.88 


8.72 


8.35 


7.93 


o.O 




0.41 


D.07 


5.76 


5.32 


— 15 




10 


1990 Aug 07 


8.83 


8.67 


8.29 


7.89 


8.0 




6.32 


5.99 


5.69 


5.24 


26 




20 


iyyu uec zz 


8.66 


8.51 


8.14 


7.71 


8.0 




6.18 


5.86 


5.54 


5.12 


— 1 




— o7 


1990 Dec 2q 


8.66 


8.50 


8.13 


7.71 


8.0 




6.17 


5.85 


5.53 


5.12 


4 




—31 


lyyl Jan zl 


8.65 


8.48 


8.12 


7.69 


7.9 




6.15 


5.83 


5.51 


5.10 


— 4 




—7 


zOOO Jul 18 










7.3 


• • ■ 


5.54 


5.15 


4.73 


4.18 









zOOO Dec Iz 










7.5 


7.67 


5.58 


5.20 


4.82 


4.29 




2 


5 


2001 Jan 17 










7.5 


7.57 


5.57 


5.21 


4.83 


4.32 







-3 


2001 Apr 12 










7.7 


7.63 


5.80 


5.45 


5.07 


4.55 




-13 


-26 


2001 May 27 










7.9 


7.90 


5.91 


5.58 


5.21 


4.71 







-14 


2001 Jun 10 










7.9 


7.78 


5.94 


5.61 


5.24 


4.75 










2001 Jun 15 










7.9 


7.71 


5.95 


5.62 


5.25 


4.76 




-4 


4 


2002 Mar 17 










7.7 


7.60 


5.93 


5.60 


5.27 


4.93 




-60 


1 


2002 Apr 30 










7.8 


7.57 


6.01 


5.68 


5.34 


4.92 




-16 


7 


2002 Jul 04 










7.8 


7.47 


5.88 


5.54 


5.20 


4.79 




-8 


2 


2002 Jul 20 










7.7 


7.46 


5.84 


5.49 


5.14 


4.73 




7 


18 


2002 Nov 04 










7.2 


7.20 


5.51 


5.12 


4.76 


4.30 




-39 


-85 


2003 Jan 11 










7.2 


7.03 


5.31 


4.89 


4.52 


4.03 







-17 



^ Photometric errors are indicated in parenthesis. References: uvby Str6mgren=LTPV project (Sterken 1983; Manfroid et al. 1991; Sterken 
et al. 1993; Spoon et al. 1994; Manfroid et al. 1995; Sterken et al. 1995); AAVSO visual=A. A. Heiden 2007, private communication; AS AS 
V-band=Pojmanski (2002); JHKL=SAAO (Whitelock et al. 1983; this work). 

^ Time delay between the epochs when spectroscopic and photometric observations were obtained (JDgpcc— JDphot)- 



18-19 between 4330-4870 A (R=2300) and 6200-7270 A 
(R=1700), and were described by Viotti et al. (1993). 

It was not possible to obtain a precise continuum nor- 
malization of the Echcllc data around Ha since the; spc;c- 
tral coverage was very small and the electron scatter- 
ing wings were very extended. Therefore, we scaled the 
Echelle spectra in flux in order to match the electron- 
scattering wings of the longslit data, as described by Lei- 
therer et al. (1994). We conservatively estimate that the 
uncertainty in the continuum after this procedure is at 
most ~ 5%. The same procedure could not be applied to 
the Echelle spectra around He I 5876 A and, as a conse- 
quence, the strength of the electron scattering wings are 
not reliable, and are likely underestimated, around this 
line. 

Additionally, the near-infrared region around 9950- 
10200 A was covered by CCD observations made on 
the 1.6m telescope of the Observatorio Pico dos Dias 
(OPD/LNA, Brazil) using the Coude spectrograph. The 
spectra typically have R ~ 10,000, and the data reduc- 
tion was accomplished using standard IRAF tasks for 
bias subtraction, flat-field correction, extraction of the 
spectrum, wavelength calibration, and continuum nor- 
malization. Those spectra were also published by Lei- 
therer et al. (1994), where further details can be found. 

2.2.2. Minimum of 2000-2003 

During this minimum we used the Fiber-fed Extended 
Range Optical Spectrograph (FEROS; Kaufer et al. 
1999) mounted on the 1.52 m telescope of ESO/La Silla 



(Chile) to obtain spectra between 2001 January and 2002 
March. The instrumental configuration provided spec- 
tral coverage between 3600- 9200 A at R=48,000, and 
the data was reduced using the FEROS pipeline (Stahl 
et al. 1999). We also used data obtained on 2003 January 
11 with the UV- visual Echelle Spectrograph (UVES) 
mounted on the ESO 8m Very Large Telescope (VLT, 
program 266.D-5655, Bagnulo et al. 2003). The resolu- 
tion provided by UVES was R = 80,000 in the wave- 
length range 3000 10400 A. 

We used the near-infrared camera CamlV on the 1.6m 
telescope of OPD/LNA (Brazil) to obtain long-slit spec- 
tra of AG Car in the wavelength range 10450 11000 A 
(i? = 7000) in 2001 June, 2002 July and 2002 Novem- 
ber. Further details on the observational setup and data 
reduction are given in Groh et al. (2007), where the data 
are also discussed. 

With the decommissioning of lUE, no ultraviolet mon- 
itoring of the star was possible during the 2000 2003 
minimum. Fortunately, a Far Ultraviolet Spectroscopic 
Explorer (FUSE) visit occurred in 2001 May 27, cover- 
ing the far-ultraviolet region between 900-1180 A [R = 
20,000). The FUSE spectrum was obtained from the 
online FUSE public-data archive. 

The spectroscopic observations obtained during mini- 
mum are summarized in Table 2. 



3. SPECTROSCOPIC VARIABILITY OF AG CAR DURING 
MINIMUM 
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TABLE 2 

Journal of spectroscopic observations of AG Car during 



MINIMUM 


Date 


Telescope 


Spectral Range (A) 


R 


1985 July 19 


lUE SWP 


1180-1900 


18,000 


1986 June 17 


1.4m ESO 


5860-5910 


60,000 


1986 June 18 


1.4m ESO 


6532-6588 


60,000 


1986 June 23 


lUE SWP 


1180-1900 


18,000 


1986 June 23 


lUE LWP 


1800-3150 


13,000 


1987 January 05 


1.4m ESO 


5860-5910 


60,000 


1987 January 05 


1.4m ESO 


6532-6588 


60,000 


1987 June 10 


1.4m ESO 


3917-3947 


60,000 


1987 June 10 


1.4m ESO 


6532-6588 


60,000 


1987 July 24 


lUE SWP 


1180-1900 


18,000 


1987 July 24 


lUE LWP 


1800-3150 


13,000 


1989 March 26 


Im CTIO/Yale 


3800-4900 


3000 


1989 December 23 


lUE SWP 


1180-1900 


18,000 


1989 December 23 


lUE LWP 


1800-3150 


13,000 


1990 April 30 


lUE SWP 


1180-1900 


400 


1990 June 16 


1.6m LNA 


9950-10200 


10,000 


1990 June 18 


1.52m ESO 


4330-4870 


2300 


1990 June 18 


1.52m ESO 


6200-7270 


1700 


1990 August 07 


lUE LWP 


1800-3150 


270 


1990 December 22 


1.6m LNA 


9970-10080 


10,000 


1990 December 28 


1.6m LNA 


6490-6710 


10,000 


1991 January 21 


1.4m ESO 


4550-4750 


60,000 


2000 July 18 


1.6m LNA 


9950-10200 


10,000 


zuuu juiy lo 




DOOU— / iUU 


iU,UUU 


2000 December 12 


1.6m LNA 


6350-7100 


10,000 


2001 January 17 


1.52m ESO 


3600-9200 


48,000 


2001 April 12 


1.52m ESO 


3600-9200 


48,000 


2001 May 27 


FUSE 


900-1180 


20000 


2001 June 10 


1.6m LNA 


10450-11000 


7000 


2001 June 15 


1.52m ESO 


3600-9200 


48,000 


2002 March 17 


1.52m ESO 


3600-9200 


48,000 


2002 April 30 


1.6m LNA 


10450-11000 


7000 


2002 July 04 


1.52m ESO 


3600-9200 


48,000 


2002 July 20 


1.6m LNA 


10450-11000 


7000 


2002 November 04 


1.6m LNA 


10450-11100 


7000 


2003 January 11 


8m ESO/VLT 


3000-10400 


80,000 



In this Section we present the spectroscopic evolution 
of AG Car during minimum epochs. Since other WN9- 
WNll stars also display variability in the spectrum on 
timescales of years (McGregor et al. 1988b; Crowther 
1997; Groh et al. 2007), investigating the variability be- 
havior of AG Car might provide insights on the evolu- 
tionary link between LBVs and WN9-WN11 stars and 
on the physical mechanism which drives the S-Dor vari- 
ability. 

3.1. Minimum of 1985-1990 

Leitherer et al. (1994) monitored the ultraviolet spec- 
trum during the minimum of 1985-1990 and concluded 
that little or no variability was present. Indeed, the low 
level of variability is consistent with the visual lightcurve 
(Fig. 1), which shows only low-amplitude photometric 
variability of ~ 0.1 — 0.3 mag during 1985-1990 (van 
Genderen et al. 1988, 1990; Leitherer et al. 1994; Stahl 
et al. 2001). The evolution of the optical spectrum from 
the end of the minimum (1990 December) until epochs 
on the rise to the lightcurve maximum (1992 June) was 
also followed by Leitherer et al. (1994). 

Here we draw our attention to the evolution of the 
optical spectrum and photometry inside the 1985-1990 
minimum. Spectra around He II 4686 A were available 
only during 1989-1991 and are shown in Fig. 2. These 
data clearly show that He II 4686 A became progressively 
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Fig. 2. — Top: AG Car spectrum around [Fe III] 4657 A and 
He II 4686 A, obtained in 1989 March (black solid line), 1990 
June (red dashed), and 1991 January (blue dotted). The resolu- 
tion of the spectra obtained in 1989 March and 1991 January were 
adjusted to the same resolution as the spectrum from 1990 June 
{R = 2300). Quoted EWs refer to He II 4686 A, and have an 
error of ~ 0.03 A. Middle: AG Car spectrum around Ha obtained 
in 1986 June (black solid hne), 1987 January (blue dotted), 1987 
June (red dashed), and 1990 December (green dot-dashed). Quoted 
EWs refer to the emission component of Ha, and have an error of 
~ 2 A wrhich is dominated by the uncertainty in the continuum 
level. Bottom: Zoom in around the variable P-Cygni absorption 
component that appears during some epochs, indicating changes 
in the ionization structure of hydrogen in the outer wind, as sug- 
gested by Crowther (1997) to explain a similar behavior seen in the 
LBV candidate He3-519. The offset of the continuum level from 
1.0 is due to the strong Ha electron scattering wings. 

weaker from 1989 until 1991, indicating a decrease of the 
effective temperature. The high-ionization line of [Fe III] 
4657 A also became slightly weaker from 1989 March 
until 1990 June. 

The Ha line was followed with a better time sampling 
than He II 4686 A (Fig. 2). It can be noticed that moder- 
ate variations were present in the emission component of 
Ha, with a decrease in the equivalent width (EW) from 
1986 June to 1987 January. The EW of Ha varied from 
-68±2 A in 1986 June to -56±2 A in 1987 January (Fig. 
2), and then remained approximately constant until 1990 
December. On the other hand, the y-band magnitude 
was 8.00 mag in 1986 June, remained roughly constant 
during 1987, and only changed significantly to 7.71 mag 
in 1990 December (Fig. 1). Therefore, the changes in 
the emission component in Ha can be explained by a net 
increase in the continuum emission compared to the line 
emission only in 1990 December - on the other hand, the 
line emission was indeed stronger in 1986 June compared 
to the rest of the minimum. 

In addition to the changes in the emission component, 
a weak P-Cygni absorption in Ha developed during some 
epochs between 1985-1990. Interestingly, a similar be- 
havior was previously noticed in the LBV candidate He 3- 
519 (Crowther 1997), whose spectrum is remarkably sim- 
ilar to AG Car's during the 1985-1990 minimum (Smith 
et al. 1994; Walborn & Fitzpatrick 2000). Crowther 
(1997) suggested that this behavior is due to changes 
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in the ionization structure of hydrogen in the outer parts 
of the wind. In this scenario, the absorption component 
seen in Ha appears when a fraction of hydrogen recom- 
bines in the outer part of the wind and disappears when 
hydrogen is mainly ionized. A similar behavior has also 
been detected in the LBV P Cygni (Pauldrach & Puis 
1990; Najarro et al. 1997). 

Comparing the evolution of the emission and absorp- 
tion components of Ha, we suggest that the stellar 
and/or wind parameters are variable even during min- 
imum epochs, and radiative transfer models are needed 
to quantify these changes. The detailed analysis of the 
1985-1990 minimum spectra will be presented in Sect. 
6.2. 

3.2. Minimum of 2000-2003 

Figure 3 displays the evolution of the optical spectrum 
of AG Car between 3690-8850 A. The variability in the 
near- infrared region from 10500 -11000 A was presented 
by Groh et al. (2007). 

In general, the spectral lines became broader from 2001 
April to 2002 March and then narrower in 2003 January, 
indicating changes in the wind terminal velocity as a 
function of temperature. That behavior is interpreted as 
due to the fast rotation of AG Car (Groh et al. 2006), but 
might also be caused by the time-dependent nature of the 
velocity field. The hydrogen Balmer and Paschcn lines 
(Fig. 3) show a significant change in strength and line 
profile inside the minimum: they became weaker from 
2001 April to 2002 March and then stronger in 2003 Jan- 
uary, indicating changes in the wind density. He II 4686 
A is very weak during the minimum of 2000 2003, with 
EW~ -100 ± 10 mA in 2001 April and 2001 June spec- 
tra, and EW~ -80±30 mA in the 2002 March spectrum. 
He II 4686 A is absent in the 2003 UVES spectrum, with 
an upper limit of 15 mA. Comparatively, this line is 
readily identified in the spectrum of AG Car during the 
1985-1990 minimum (Fig. 2; Stahl 1986; Viotti et al. 
1993; Leitherer et al. 1994; Smith et al. 1994; Walborn 
& Fitzpatrick 2000: Stahl et al. 2001). 

Meanwhile, the He I, Fc III, and N II lines became 
weaker from 2001 to 2003, while the very weak Fc II 
and Mg II lines seen in the 2001 and 2002 spectra are 
prominent in the data from 2003. This behavior indicates 
that the effective temperature was decreasing from 2001 
to 2003. Therefore, the above picture is consistent with 
the star moving out of minimum from 2001 to 2003, as 
suggested by the visual lightcurvc (Fig. 1). 

Interestingly, multiple absorption components ap- 
peared during 2001-2003 in the hydrogen, Fe II, and Ni II 
lines (Fig. 4), and are reminiscent of what was seen dur- 
ing 1991 1993 (Leitherer et al. 1994; Stahl et al. 2001) 
when the star was also moving towards maximum. Three 
absorption components in the Fc II lines can be securely 
identified in the data from 2000 2003, and they are cen- 
tered at -72, -101, and -148 kms^^ in the FEROS 
spectra from 2001-2002. These absorption components 
change in velocity and in strength in timescales as short 
as 1 year, as first noted by Stahl et al. (2001). Our 
data also present a similar behavior, and in the 2003 
UVES spectrum, these absorption components shift to 
—74, —112, and —153 kms~^. An additional absorption 
component centered at —202 kms~^ is also present in 



the Fe II and hydrogen Balmer lines in 2003. 

The —72 kms^^ component could arguably be formed 
in the distant circumstellar nebula, which has an expan- 
sion velocity of ~ 70 kms~^ (Nota et al. 1992), but that 
component is not observed at all epochs (Leitherer et al. 
1994; Stahl et al. 2001, this work). The other multi- 
ple absorption components are unlikely to be formed in 
the distant circumstellar nebula of AG Car, since they 
change in velocity and are approximately centered at the 
wind terminal velocity found in previous epochs. Thus, 
the multiple absorption components in AG Car are likely 
not reminiscent of the multiple ejecta absorption compo- 
nents seen in the massive LBV Eta Carinac (Gull et al. 
2005; Nielsen et al. 2005; Johansson et al. 2005; Gull 
et al. 2006). We support the scenario, proposed by Stahl 
et al. (2001), that the multiple-absorption features are 
formed much closer to the star and are due to regions 
in the wind which have a density decrease (or enhance- 
ment) due to the temporal variability of the star. Indeed, 
this behavior is predicted by detailed radiative transfer 
models which include time-dependent effects (Groh et al. 
2008; J. H. Groh et al. 2009, in preparation). 

4. THE MODEL 

To analyze the spectra of AG Car, we used the radia- 
tive transfer code CMFGEN (Hilher 1987, 1990: HiUier 
& Miller 1998, 1999; Busche & HiUier 2005), which has 
been successfully used to model the spectrum of LBVs 
and related objects (Smith et al. 1994; Najarro et al. 
1994, 1997, 2009; Najarro 2001; Hillier et al. 1998, 2001, 
2006; Figer et al. 1998; Drissen et al. 2001; Bresolin et al. 
2002; Marcolino et al. 2007). The code has been exten- 
sively discussed in the aforementioned references, and be- 
low we concisely describe the main characteristics of the 
code. 

CMFGEN assumes an spherical-symmetric, steady- 
state outflow, and computes continuum and line forma- 
tion in the non-LTE regime. Each model is defined by 
the hydrostatic stellar radius i?*, the luminosity L^,, the 
mass-loss rate M, the wind terminal velocity v^o , the stel- 
lar mass M, and by the abundances of the included 
species. We define the hydrostatic radius R^, as the ra- 
dius where the velocity is i;(i?*) = fsonic/3, in order to 
avoid any effects due to the strong wind in determination 
of i?*, which is the case for v > Wsonic/3. Thus, the veloc- 
ity of the hydrostatic core, where ii* and are defined, 
is typically 4 — 5 kms~^ for AG Car during minimum. 

CMFGEN docs not solve the hydrodynamic equations 
of the wind, and therefore, a velocity law v{r) has to 
be assumed a priori. In CMFGEN, the velocity law is 
parameterized either by a single beta-type law or by a 
two-beta velocity law which is modified at depth to reach 
a hydrostatic structure. The stellar mass M is needed to 
compute the density structure below and close to the 
sonic point, and therefore changes in M might affect the 
spectroscopic analysis. In this paper we assume M = 
70 Mq, since lower values of M would imply that the 
star is above the Eddington limit modified by rotation 
(Paper II). A thorough discussion on the impact of M 
on the spectral analysis will be provided in Paper II. 
The velocity structure below the sonic point is iterated 
to fulfill the wind momentum equation within 10% for 
t; < 8 kms~^. A beta- type law is joined to this structure 
just below the sonic point. Although the code solves for 
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Fig. 3. — Evolution of the optical spectrum of AG Car in the region 3690-8850 A during the 2000-2003 minimum (see text for discussion) . 
The spectra were obtained in 2001 April (black solid line), 2002 March (blue dashed line), and 2003 January (red dot-dashed line). Data 
obtained from 2000 July until 2001 Juno show no significant variability. The spectrum from 2003 January has a small gap in wravelength 
coverage between 8540 — 8660 A, while Ha is saturated and not shown. 



a hydrostatic structure, there might stiU be uncertainties 
in the density structure below the sonic point due to the 
proximity to the Eddington hmit and the choice of M, 
rotation rate, and viewing angle. 

CMFGEN allows for the presence of clumping within 
the wind using a volume-filling factor approach. The 
wind is assumed to be homogenous close to the star, be- 
come partially clumpy at a characteristic velocity scale 
and achieves full dumpiness, with a volume-filling 
factor /, at large distances, as follows: 



/(r) = / + (l-/)cxp[-i;(r)K 



(1) 



Full line blanketing is included consistently in CMF- 
GEN through the concept of super levels, which groups 
similar energy levels into one single superlevel to be ac- 
counted for in the statistical equilibrium equations. In 
the analysis of AG Car, our initial atomic model con- 
sisted of H, He, C, N, O, and Fe in order to constrain the 
parameter regime to be explored. For the final modeling 
we included additional elements and ionization stages. 
The atomic models adopted for the final modeling are 
summarized in Table 3, and included H, He, C, N, O, 



Na, Mg, Al, Si, Cr, Mn, Fe, Co, and Ni. 

The synthetic spectra of AG Car presented in this pa- 
per were computed using CMF_FLUX (HiUier & Miller 
1998; Busche & HiUier 2005) and did not include effects 
due to rotation. As discussed in Busche & HiUier (2005), 
rotation mainly affects spectral lines formed close to the 
photosphere. Specifically in the case of AG Car during 
minimum, only Si IV and He II lines have their profile 
significantly changed by rotation (Groh et al. 2006; Pa- 
per II). 



5. DISTANCE AND SYSTEMIC VELOCITY 

The spectroscopic analysis using CMFGEN cannot 
constrain the distance to AG Car. Throughout this 
work, we assume that AG Car is located at c? = 6 kpc 
(Humphreys et al. 1989; Hoekzema et al. 1992). 

Stahl et al. (2001) re-determined the heliocentric sys- 



temic velocity of AG Car and found Vs- 



10±5 kms' 



suggesting a lower kinematical distance in the range of 
5-6 kpc, instead of a kinematical distance of 6.0-6.5 kpc 
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Fig. 4. — Evolution of the multiple absorption components seen 
in the AG Car spectrum during the minimum of 2000-2003. To 
illustrate their general behavior, we display the Fe II 5170 A line 
(top) and H7 4340 A (bottom). The high- velocity multiple ab- 
sorption components are likely caused by time-dependent effects, 
while the —72 kms~^ component may also be formed in the cir- 
cumstellar nebula of AG Car, which has an expansion velocity of 
~ 70 kms-l (Nota et al. 1992). 

from Humphreys et al. 1989^. Taking the errors into ac- 
count, both determinations are consistent (Stahl et al. 
2001). We also measured the centroid of the [Fe II] and 
[N II] lines in our high-resolution data and obtained Vgys 
similar to the value found by Stahl et al. (2001). 

A note of caution regarding the kinematical distance 
of early- type stars should be pointed out: the distance 
to young stellar clusters based on spectroscopic parallax 
are systematically lower than their kinematical distance 
(e.g., Blum et al. 2000, 2001; Figueredo et al. 2002, 2008). 
In addition, about 10-30 % of Galactic 0-type stars are 
runaways (Gies 1987; de Wit et al. 2005), implying that 
their derived kinematical distance might be meaningless. 
However, in the case of AG Car, the empirical luminosity- 
amplitude relationship found for LBVs (Wolf 1989) and 
the distance versus reddening relationship of nearby stars 
both support a distance of 6 kpc (Humphreys et al. 1989). 

The bolometric luminosity L^,, stellar radius -R*, and 
mass-loss rate M obtained through the spectroscopic 
analysis are distance-dependent, while the effective tem- 
perature and chemical abundances are only weakly, if at 
all, dependent on d. Following Schmutz et al. (1989), 
Najarro et al. (1997), and Hillier et al. (1998), the phys- 
ical parameters of AG Car obtained in this work (see 
Sect. 6) should be scaled to other distances according 
to: L^, (X cP, M (X d^-^, iZ* oc d, and Tcff oc d° . 

6. RESULTS 

6.1. Surface chemical abundances 

The surface chemical abundances obtained for AG Car 
are listed in Table 4. The values obtained are consistent 
when comparing different epochs from both minima, sup- 
porting the physical parameters derived for those epochs 

^ We updated the values of the kinematical distance presented by 
Humphreys et al. (1989), rf^in ~ 6.4 — 6.9 kpc, in order to account 
for a revised galactocentric distance of the Sun of Ro = 8 kpc and 
©0 = 220 kms"'^, obtaining dj^in ~ 6.0 — 6.5 kpc, depending on 
whether the Galactic rotation curve of Fich et al. (1989) or Brand 
& Blitz (1993) was chosen. 



TABLE 3 
Final atomic model used 
in the analysis of 
AG Car 



on 








HI 


20 


30 


435 


He I 


40 


45 


233 


He II 


22 


30 


435 


CII 


16 


54 


366 


cm 


30 


54 


268 


CIV 


21 


64 


1446 


NI 


44 


104 


885 


Nil 


41 


84 


530 


NIII 


30 


63 


344 


OI 


69 


161 


358 


on 


19 


57 


372 


OIII 


17 


79 


450 


Nal 


22 


71 


1550 


Mg II 


18 


45 


358 


Al II 


21 


65 


1348 


Si II 


22 


43 


182 


Si III 


20 


34 


92 


Si IV 


22 


33 


183 


Cr II 


42 


310 


2127 


Cr III 


40 


209 


4077 


Cr IV 


29 


234 


4466 


Cr V 


30 


223 


3207 


Mn II 


33 


333 


1348 


Mn III 


33 


333 


4038 


Mn IV 


39 


464 


11103 


Mn V 


16 


80 


676 


Fell 


61 


261 


3837 


Felll 


33 


191 


6935 


FelV 


58 


383 


2816 


Fe V 


25 


168 


1938 


Co II 


54 


421 


3267 


Co III 


22 


174 


3837 


Nill 


29 


204 


4728 


Ni III 


28 


220 


5181 


Ni IV 


36 


200 


2337 


Ni V 


46 


183 


1524 


Note. 


— For each 


species, 



an additional one-level ion 
corresponding to the high- 
est ionization stage was in- 
cluded explicitly in the rate 
equations (i.e. C V for car- 
bon, for instance), but are 
not shown above for brevity. 
The columns correspond to 
the number of included su- 
per levels (Afg), total energy 
levels (Afp), and the number 
of bound-bound transitions 
included for each ion (A^x). 

(Sect. 6.2 and 6.3). 

The helium abundance of AG Car was determined as 
He/H=0.43 ± 0.08 (in number) in order to reproduce the 
intensity of He I, He II and hydrogen lines in the opti- 
cal and near infrared spectrum. This value is in agree- 
ment with previous works (Smith et al. 1994; Leitherer 
et al. 1994; Stahl et al. 2001), and with the helium con- 
tent for the LBV phase predicted by evolutionary models 
(Meynet & Maeder 2000). Unhke Eta Carinae (HiUier 
et al. 2001) and HD 316285 (Hillier et al. 1998), the de- 
duced helium abundance of AC Car does not scale with 
the derived mass-loss rate. This is because AG Car simul- 
taneously presents He I and He II lines in the spectrum 
during minimum, while Eta Car and HD 316285 do not. 
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In addition to the helium content, it is essential to de- 
termine the abundanee of the CNO-elements to eonstrain 
the evolutionary stage of AG Car. A nitrogen mass frac- 
tion of (7.2 ± 2) X 10^3 (11.5 ± 3.4 times the solar value; 
hereafter, we used the solar abundance values listed in 
Grevesse et al. 2007 and references therein) was obtained 
based on numerous N II lines in the optical spectrum 
(4601-4643, 6482 A) and on ultraviolet N III lines. The 
carbon abundance was determined based on the strength 
of optical C II Hues at 6578-6582 A and 7231-36 A, yield- 
ing a carbon mass fraction of (2.4±0.7) x 10"'' (0.11±0.03 
times the solar value). The errors in the carbon and ni- 
trogen abundance are about 30%. The oxygen mass frac- 
tion was determined to be (2.4 ± 1) x lo'"" (0.04 ± 0.02 
times the solar value) using the O I lines at 7774 8446 
A and O II 4620 A. As it was not possible to fit both 
O I lines with the same model, the uncertainty in the 
oxygen abundance is around 50%. Nonetheless, it is safe 
to point out that the oxygen abundance on the surface 
is very low compared to the solar value. Together with 
the high nitrogen content and carbon depletion, we con- 
clude that there is certainly CNO-processcd material on 
the surface of AG Car. The total C-l-N-l-0 abundance of 
AG Car is in broad agreement with a solar composition. 

The CNO cycle is also supposed to enhance the Na 
abundance in massive stars by a factor of ~ 4 (Prantzos 
et al. 1986), which could be investigated using the wind 
component of the optical resonance lines of Na I 5885- 
5890 A. We found that Na is overabundant by a factor 
of 3.3 ± 1 in comparison to the solar value. 

The rich emission-line spectra of AG Car also allowed 
us to constrain the Si and Mg surface abundances. The 
former was obtained using the ultraviolet and optical 
lines of Si II, Si III, and Si IV, while the latter was con- 
strained using Mg II lines at 4481, 10915, and 10952 
A. For both species, the derived abundance is consistent 
with the solar value within 30%. 

Iron contributes significantly to the emergent spectrum 
of AG Car and, together with other iron-group elements, 
is the main responsible for the severe line blanketing seen 
in the ultraviolet spectrum. Solar abundance was ini- 
tially assumed for Fe, Ni, Cr, Mn, and Co, and was found 
to be consistent with the line blanketing seen in the UV. 
The severe line blending with Fe lines and uncertainties 
in the atomic parameters make it impossible to constrain 
the Ni, Cr, Mn, and Co abundance better than a factor 
of 2. The Fe content can be better constrained to within 
±30% of the solar value using the numerous Fe II and 
[Fe II] lines in the optical spectrum and the Fe II forest 
in the ultraviolet. 

6.2. Stellar and wind parameters during the 1985-1990 

minimum 

We divided the data obtained from 1985 1990 into 
three siibsets which were analyzed separately as follows. 
The first epoch comprised spectra from 1985 July un- 
til 1986 June, since the Ha spectrum obtained in 1986 
June showed significantly enhanced emission compared 
to other epochs (Fig. 2). In the second subset, the data 
from 1987 January until 1990 June were averaged and 
analyzed as a single-epoch observation, since very little 
variability was present. We believe that the fundamen- 
tal parameters obtained during these epochs are the least 



TABLE 4 

Surface chemical abundances of AG Car 



Species Number fraction Mass fraction Z/Zq'^ 

H 
He 
C 
N 
O 
Na 
Mg 
Al 
Si 
Cr 
Mn 
Fc 
Co 
Ni 



^ Ratio between the mass fraction of a given species in AG Car 
and the respective solar mass fraction. Quoted errors indicate 
the range of abundances consistent with the observations. Solar 

abundances arc from Grcvcssc ct al. (2007). 

affected by time-dependent effects. The third subset con- 
sisted of data from 1990 December until 1991 January, 

since the He II 4686 A emission was noticeably weaker 
than other epochs during this minimum (Fig. 2). 

Figures 5, 6, and 7 show the comparison between the 
best CMFGEN models obtained for the 1985-1990 mini- 
mum phase and the ultraviolet, optical, and near-infrared 
spectra of AG Car. The models quantitatively reproduce 
the observed spectrum with a superb fit in most spectral 
regions, in particular in the ultraviolet, where the line 
blending and blanketing are severe. 

Table 5 presents the fundamental parameters of 
AG Car during the 1985-1990 minimum, and below we 
discuss the diagnostics used to derive each of them. 

6.2.1. Effective temperature (Tes) and temperature at the 
base of the wind (Ti,) 

During this minimum, the most sensitive diagnostic for 
deriving T^s was the ratio between He II 4686 A and He I 
hues at 4471 A, 4713 A, 5876 A, 6678 A, and 10024 A. 
Changes of a few hundred K in Teff were sufficient to 
change the EW of He II 4686 A by a factor of 2. In ad- 
dition, the ionization structure of C, N, and Si were also 
sensitive to Tcff and yielded similar values to Toff. For 
the C ionization structure, we used the ultraviolet C II 
1331-1334 A, C III 1175 A, and C III 1247 A and the op- 
tical C II 6578-6582 A and C II 7231-7236 A diagnostic 
lines. In the case of the N ionization structure, we com- 
pared the strength of several optical N II lines such as 
N II 4601-07-30-43 A, N 11 6482 A, and N II 6613 A with 
the N III triplet at 1747-1749-1751 A. The absence of 
N III Hues aroimd 4500 A was used to constrain the max- 
imum value of Toff during 1987-1990. The Si ionization 
structure was constrained through the relative strength 
of Si II, Si III, and Si IV lines, using lines of Si II 6347 
A, Si III 1280-1290 A, Si III 4553-4568-4575 A, Si IV 
1393-1402 A, and Si IV 4088-4116 A. 

Using such a large number of relatively high-sensitive 
diagnostics available for T^s, we obtained ~ 
22, 800± 500 K during 1987-1990. There was no evidence 
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from the 1985 July and 1986 June UV spectra that Tcff 
was higher than during 1987-1990 within the errors, and 
therefore we also adopted T^s ~ 22, 800 ±500 K for both 
epochs. Even if the He II 4686 A emission was slightly 
weaker in 1990 June than in 1989 March (Fig. 2), chang- 
ing Teff by a few hundreds of K within the given error is 
sufficient to fit this line. Thus, for simplicity and due to 
the absence of other extremely sensitive diagnostic lines, 
we assume that T^s was constant during 1987-1990. The 
smaller He II 4686 A emission during 1990 December- 
1991 January, when AG Car was brightening, suggests 
that the effective temperature was approximately 1000 K 
lower at that epoch {T^s ~ 21, 500 ± 500 K). 

Since the photosphere of AG Car is extended, we found 
that Teff depends on the derived value of Af (Sect. 6.2.2), 
wind velocity law (Sect. 6.2.3), and M (^ 70 Mq, Paper 
II). In addition, we found that the velocity and density 
structure below the sonic point significantly affect the 
value of Teff- In our best models, the velocity structure 
was iterated at depth in order to fulfill the momentum 
equation below the sonic point. Initial models, which 
assumed a fixed scale height across the extended atmo- 
sphere, yielded higher values of Tcff by up to ~ 1000 
K. 

The value of is defined in this paper as the tem- 
perature of the hydrostatic core (i.e at the base of the 

wind, Ri,) where v = i'sonic/3 (Sect. 4), and the differ- 
ence between T-^ and Tdi of a given epoch is directly 
related to the cxt(^nsion of the atmosphere. We ob- 
tained T*=26,200 ± 500 K during 1987 January - 1990 
June and ~ 24, 640 ± 500 K during 1990 December- 
1991 January. Since we assumed that Tgff was constant 
from 1985 July-1990 June, the higher mass-loss rate of 
AG Car obtained from 1985 July-1986 June compared 
to 1987 January 1990 June (see Sect. 6.2.2) implies 
that the atmosphere was more extended between 1985 
July-1986 June. Therefore, T^, was slightly higher be- 
tween 1985 July-1986 June (T* ~ 26, 450 ± 500 K) than 
during 1987 January - 1990 June. 

6.2.2. Mass-loss rate and wind clumping 

The mass-loss rate was obtained by reproducing the 
intensity of the strongest lines in the observed spec- 
trum of AG Car: namely, those of hydrogen, helium, 
and nitrogen. Using those diagnostics we obtained M = 
1.5 X 10~^ Mqyt-^ during 1987-1990. Spectra from 1990 
December-1991 January are also consistent with the 
same M. Figure 8 illustrates the effects of wind clump- 
ing on selected spectral lines of AG Car, presenting CMF- 
GEN models assuming different volume-filling factors but 
with the same value of M/""-^ = 4.7 x lO"'"' Mgyr-^and, 
thus, equally reproducing the amount of Ha emission. 
We derived a volume-filling factor of / = 0.1 based on 
the strength of the electron-scattering wings of Ha, H/3, 
H7, H(5, and He I 6678 A. 

Similar to what has been found for 0-stars (Hillier 
et al. 2003; Bouret et al. 2003, 2005), we found that 
clumps must be created very deep in the wind, close to 
the sonic point (i.e. Vc ~ 20 kms~^), in order to repro- 
duce the weak electron-scattering wings of He II 4686 A. 
The presence of clumps very close to the stellar surface 
is in line with the latest polarimetric results by Davies 
et al. (2005, 2006, 2007, 2008). 



The statistical error in M when using a single spec- 
tral line, such as Ha, as diagnostic is as small as 5%. 
However, taking into account that multiple lines need to 
be fitted, that additional parameters can also change the 
line strength and have their own errors, and the assump- 
tions in the model (prescribed velocity law, steady-state, 
spherical-symmetric wind, etc.), we estimate that the un- 
certainty in M amounts to 30%. We found that models 
assuming either / < 0.05 or / > 0.25 do not provide 
reasonable fits to the electron-scattering wings present 
in the AG Car spectrum (Fig. 8). 

As shown in Sect. 3.1, moderate spectroscopic vari- 
ability is seen during the minimum only in the Ha profile 
obtained in 1986 June, which was ^ 20% stronger than 
during the rest of the minimum (1987-1990). To inves- 
tigate the nature of the changes seen in Ha and explain 
the behavior of AG Car during 1985-1990, we explored 
the parameter space around the best CMFGEN radia- 
tive transfer model obtained for 1987-1990 (see Table 
5). This model spectrum does not have a P-Cygni ab- 
sorption component in Ha and provides a nice fit to the 
Ha profile from 1987 June - 1990 June. 

We found that the Ha emission in 1986 June can be 
explained by an increase of ~ 30% in M compared to 
the value obtained for the rest of the minimum (M = 
1.5 X 10~^ M0yr^^, see Sect. 6.2.2), corresponding to 
M = 1.9 X lO"-"' MqYT-^. However, increasing M by 30% 
is not sufficient to change the H ionization structure and 
does not produce a weak P-Cygni absorption component. 
Furthermore, the change in Af by 30% in the model pro- 
duces a modest 1% increase in the y-band flux, or 0.01 
mag. Interestingly, such an amount of photometric vari- 
ability induced by the change in M is compatible with 
the microvariability detected in AG Car during minimum 
(van Genderen et al. 1988, 1990). 

As has been suggested in order to explain the spectro- 
scopic changes seen in the LBVs He 3-519 (Crowther 
1997) and P Cygni (Pauldrach & Puis 1990; Najarro 
et al. 1997), we also found that the H ionization struc- 
ture of AG Car could be affected by changes in M, Ri,, 
and Li,. However, in the parameter regime found during 
1985-1990, we determined that even if M is increased by 
a relatively large factor of 4, the model spectrum does 
not show a weak P-Cygni absorption component, as seen 
in the observations. Therefore, a decrease in Teff and/or 
is required in order to fit the weak P-Cygni absorp- 
tion in the 1986 June spectrum. The decrease in T^s is 
not supported by the ultraviolet spectrum taken in 1985- 
1990, while the change in would produce, according 
to our models, significant photometric changes of at least 
'~ 0.5 mag, which is not observed during 1985 1990. 

6.2.3. Wind terminal velocity and velocity law 

We used the maximum velocity of the saturated P- 
Cygni absorption component of the UV resonance lincs^*^ 
of C II 1334-1335 A, Si IV 1393-1402 A, and Mg II 2798- 
2810 A as a diagnostic to derive the wind terminal ve- 
locity. In the case of C II 1334 1335 A both lines arc 
blended, and only the blue component was used. The 
determination of Wqo from UV resonance lines is very 

The so-called Ublacki see Prinja et al. (1990). 
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Fig. 5. — Comparison between the ultraviolet spectrum of 
AG Car from 1989 December 23 (black line) and the best CM- 
FGEN model for that epoch reddened using E{B — V)= 0.65 and 
Rv = 3.5 (red hne) in the spectral region between 1260 — 3150 A. 
The same model provides similarly good fits to the other ultra- 
violet spectra of AG Car taken between 1985-1990, since there 
was little variability in the UV spectrum of AG Car during these 
epochs (Leitherer et al. 1994). Line identification for the strongest 
spectral features is presented. The spectrum in the region between 
1500 — 1950 A is highly blanketed and dominated by blends of 
Fe III, Fe IV, and Ni IV lines. Notice that Fe II lines are absent. 
The P III doublet lines at 1344.32-1344.85 A were not included in 
the models. 

weakly dependent on the assumed velocity law. Since the 
UV resonance lines essentially have no variability dur- 
ing minimum (Leitherer et al. 1994), all epochs between 
1985-1989 yield = 300 ± 30 kms-^. 

Additional constraints on the value of Uoo can be ob- 
tained from the FWHM of the emission component of 
the strongest optical lines, such as Ha, H/3, H7, and 
He I lines, and from the P-Cygni absorption profile of 
the He I lines. Using these diagnostics, we derived 
Voo — 270 ± 50 kms~^, which is slightly lower but still in 
line with the value obtained from the UV resonance lines. 
The optical lines provided only a lower limit to Voc , since 
they were much more sensitive to the adopted velocity 
law and to the density structure of the wind than the 
UV resonance lines, which explains the different results 
comparing optical and UV analyses. Using only the op- 
tical lines, lower values of v^o were obtained when low 
values of the wind acceleration parameter f3 were used. 
In particular, a normal velocity law with /3 = 1 would 
require Vqo as low as 230 kms~^, while higher values 
of /3 such as 3-4 would require a value of v^o closer to 
that obtained from the UV lines. Thus, ultimately, we 
favor the use of the value obtained from the UV lines, 
Voo = 300 ± 30 kms~^, and that (3 is roughly 3. A high 
value of P is supported by several other diagnostics. Low 
values of beta did not simultaneously provide reasonable 
fits to UV resonance lines and to the optical H I and He I 
line profiles, in particular to their absorption profile. In 
addition, the ratio between Ha and the other Balmer 
lines such as H/3 and H7 required /3 > 3. Values of /3 < 3 
also yielded too much radiative acceleration beyond the 
sonic point (Paper II). 

Evidence of an extended acceleration zone in the outer 
wind is suggested by the wind hydrodynamics: a velocity 



law using a single value of f3 yielded too much radiative 
acceleration for distances greater than ~ 10^ R^^ (Paper 
II). However, we chose to use a usual beta- type law in 
order to not introduce extra free parameters in the anal- 
ysis. 

Line-driven winds are intrinsically unstable, and tur- 
bulent motions are expected to arise (Owocki et al. 1988; 
Feldmeier 1995). Such a turbulent velocity field pro- 
duces a non-saturated absorption wing extending to ve- 
locities higher than Voo on most UV resonance lines in 
the AG Car spectrum. In addition, a redshift of the 
emission components of the stronger H I, He I, and Fe II 
lines is also produced due to the turbulent velocity field 
(HiUier 1987, 1989; Catala et al. 1984). In CMFGEN, it 
is possible to account for a microturbulent velocity Wturb, 
which was set to 20 kms^^. The use of higher values 
of Yturb causes a redshift of the emission component of 
the stronger He I and hydrogen lines, which is not seen 
in the observed spectrum. This also explains the shift 
of ~ -|-10kms~^ of the Ha emission component, rela- 
tive to the systemic velocity, measured by Stahl et al. 
(2001) (the same behavior is observed in our data). We 
obtained satisfactory fits to the extended, non-saturated 
absorption wing using vturb = 20 kms^^. 

Leitherer et al. (1994) noticed that during minimum 
the UV resonance lines of Si IV 1393-1402 A have an 
extended, non-saturated absorption wing with velocities 
up to ~ —700 kms^^. Those authors interpreted this 
component either as evidence of a fast polar outflow of 
AG Car or as being caused due to blending with Fe 
lines (following Hubeny et al. 1985). Our models pre- 
dict severe blending due to Fe III and Fe IV around both 
Si IV lines, which can mimic a high- velocity absorption 
wing. Therefore, we do not support the idea that the 
blueshifted, high-velocity component seen in the Si IV 
lines are actually due to a fast polar wind; instead, we 
support that they are caused by Fe III and Fe IV absorp- 
tion lines. The high projected rotational velocity inferred 
for AG Car by Groh et al. (2006), and the inference that 
we are seeing AG Car close to edge-on, provides further 
evidence that the high- velocity absorption does not come 
from the polar wind. 

6.3. Stellar and wind parameters during the 2000-2003 

minimum 

Unlike the previous minimum, the observations of 
AG Car between 2000-2003 show significant temporal 
variability. The data were divided in four datasets 
consisting of data between 2000 July-2001 June, 2002 
March-2002 July, 2002 November, and 2003 January. 
There was little variability during the time intervals of 
2000 July-2001 June and 2002 March-2002 July; there- 
fore, the average spectrum was used for those inter- 
vals. As a consequence of the larger spectral coverage 
and higher spectral resolution, it was possible to use a 
larger number of optical diagnostic lines during 2000- 
2003 in comparison to the previous minimum. The spec- 
tra taken before 2000 July and after 2003 January indi- 
cate Teff < 13,000K, constraining the duration of this 
minimum to at most 2.5 years. The maximum T^s was 
only achieved during ^^1-1. 5 years (Table 5). 

The absence of ultraviolet observations did not affect 
the determination of the stellar and wind parameters of 
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Fig. 6. — Optical spectra of AG Car taken at different epochs 
during the 1985-1990 minimum (black line) compared with the 
respective best CMFGEN model for each epoch (red line; see Table 
5). Note the presence of [Fe 111] 4658 A during 1989 March in panel 
h. 



AG Car since a much larger spectral coverage in the opti- 
cal and near-IR was available compared to the 1985-1990 
minimum, providing enough diagnostics. In addition, the 
main parameter obtained from ultraviolet lines, namely 
the wind terminal velocity, was determined based on the 
numerous P-Cygni absorption components of He 1, N 11, 
Fe II, and hydrogen Balmer and Paschen lines. For the 
data taken in 2001 May, it was also possible to constrain 
Woo using C 111 1173 A from the FUSE spectrum, yielding 
a similar value of Voo as the one found from the optical 
lines. 

The optical diagnostic lines used to obtain the stellar 
and wind parameters are similar to those described in 
Sect. 6.2 and are just outlined here. The value of M 
was obtained in order to reproduce the strength of the 
H, He, N, O, Mg, Si and Fe lines, while the volume- 
filling factor / was constrained using the intensity of the 
electron-scattering wings of the strongest H, He 1, and 
Fe II lines. 

The value of was determined based on the ion- 
ization structure of C, N, O, and Si, and on the rela- 
tive strength between H and He I lines, assuming the 
He abundance obtained for the previous minimum. The 
presence of very weak He II 4686 A emission was an ad- 
ditional constraint for the maximum value of Tcs. The 
N ionization structure was constrained by comparing the 
strength of numerous spectral lines^^ of N I (e.g., 7442, 

We refer to the atlas of Stahl et al. (1993) for a comprehensive 
identification of the spectral lines found in the parameter regime 
investigated in the present paper. 
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Fig. 7. — Optical and near-IR spectra of AG Car taken at differ- 
ent epochs during the 1985-1990 minimum (black line) compared 
with the respective best CMFGEN model for each epoch (red line; 
see Table 5). Notice in panel g that He II 4686 A is broader in 
the observations than in the CMFGEN model, which is due to the 
fast rotation of AG Car (Groh et al. 2006; Paper II). In addition, 
the observed He I 4713 A line profile is broader than the model in 
1991 January (panel g). 

7468, 8184-88, 8680-83-86, 8703-11-18 A), N II (e.g., 
4607-4630-4643, 6482 A), and N 111 4097 A. Likewise, 
several O I (e.g. 8446 A) and many weak O 11 lines (e.g., 
3973, 4349, 4070-72, 4317-19-49, 4649 A) present in the 
optical spectrum of AG Car were used to derive Toff- As- 
suming the C abundance obtained for the previous min- 
imum, the C ionization structure could be constrained 
based on the strength of several C II lines, such as those 
at 6578 - 82 A and 7231 - 36 - 37 A. Optical spectral 
hues of Si 11, Si 111, and Si IV were also used as diag- 
nostics of Toff- The values of Tes provided by different 
diagnostics were consistent within ^ 500 K. 

The strength of Fe III and some Fe 11 lines were well re- 
produced by the CMFGEN models, supporting the value 
obtained for T^s and M. Note that a few Fe II lines, in 
particular those formed in the outer part of the wind such 
as Fe II 5169 A, were stronger in the models than in the 
observations. We suspect that these lines are affected by 
time-dependent effects. 

Figures 9a, 9b, and 9c present the observed AG Car 
spectra from 2001 April, 2002 March, and 2003 Jan- 
uary, respectively, and the corresponding best CMFGEN 
model for each epoch. In general, a reasonably good fit 
to the observed spectrum could be obtained. However, in 
the 2001 April data a weak excess emission on top of the 
P-Cygni profile of H and He I is present approximately 
between -1-130 and -1-250 kms""'^, and is not reproduced 
by the models. Such an excess emission could be caused 
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Fig. 8. — Optical spectra of AG Car taken at different epochs 
during the 1985—1990 minimum (black line) compared with CMF- 
GEN models assuming different volume-filling factors / and with 
M scaled to produce a similar amount of Ha emission (i.e., all 
models have the same Mf"^-^ = 4.7 X 10~^ MQyr~^). For each 
epoch we display models with / = 0.05 (dotted blue line), / = 0.10 
(dashed red line), and / = 0.25 (dot-dashed green line). Notice 
that the model with / = 0.25 produces too strong electron scat- 
tering wings and increased He II 4686 A emission. Increasing the 
characteristic velocity at which clumps start to form (vc) from 20 
to 70 kms~^ produces a similar effect to the He II 4686 A emission 
as increasing the volume-filling factor, since this line is formed very 
deep in the wind. Note in panel c that the Si IV 4088 - 4116 A 
absorptions are broader in the observations than in the CMFGEN 
model because rotational broadening was not included (Groh et al. 
2006; Paper II). 

by a latitude-dependent wind or by a velocity law differ- 
ent than the usual beta-type law assumed in this work. 
The values obtained for T^,, Tcs, Voo, M, and / are listed 
in Table 5. 

6.4. The FUSE spectrum of AG Car and the extreme 
line blanketing in the far ultraviolet 

Since the launching of the Far Ultraviolet Spectro- 
scopic Explorer (FUSE) satellite, high-resolution far- 
ultraviolet spectroscopy of OB and WR stars has pro- 
vided invaluable constraints on their fundamental prop- 
erties, since a significant fraction of the stellar light is 
emitted in this spectral region and important diagnostic 
lines fall in the far ultraviolet. For a recent review on 
the accomplishments of the far-UV research of massive 
stars, see Crowther (2006). 

Figure 10 shows the comparison between the archival 
FUSE spectrum of AG Car obtained during 2001 May 
and the best CMFGEN model for that epoch reddened 
using E{B -V) = 0.65 and Ry = 3.5 (see Sect. 6.5). 
The CMFGEN model spectrum matches reasonably well 
the observations, and the FUSE spectrum is consistent 
with the presence of a dense wind, indicating M in the 
range 3 — 4 x 10~^ M0yr~^. The relative strength of 
Fe II and Fe III lines is the main diagnostic for T^g, and 
agrees with the value derived from the optical/near-IR 
analysis. The many Fe II and Fe III transitions provide 
a good constraint on Voo, as does C III 1173 A: all these 
diagnostics indicate a very low wind terminal velocity 
during 2001 May {v^ ~ 90 - 130 kms'^). 

It can readily be noticed in Fig. 10 that line blanket- 
ing is extreme in the far-UV, and the determination of 
the continuum level with confidence is impossible from 



the observations alone. Besides C III 1173 A, the iden- 
tification of other spectral features is hampered by the 
severe line overlapping present in the far-UV. In order to 
determine which species are responsible for a given spec- 
tral feature, we computed the observed spectrum due to 
a given ion using CMF.FLUX (Busche & Hillier 2005). 
We found that most of the observed spectral features in 
the far-UV are actually a blend of line transitions, even 
in the case of C III 1173 A, which is contaminated by 
Fe II and Fe III hues (Fig. 10). We present in Fig. 10 
the separate contributions of lines of a given ion to the 
observed spectrum of AG Car. Most of the line blanket- 
ing in the FUSE range is provided by Fe III, while the 
contribution of other species are localized to a certain 
spectral range. In the 2001 spectrum of AG Car, most 
of the Fe II blanketing in the far-UV is in the spectral 
region 1050-1180 A; Cr III contributes at 1000-1070 A, 
Co III at 920-965 A, and Ni III around 950-980 A. 

As has been seen in Galactic OB-type stars (Pellerin 
et al. 2002), the interstellar absorption due to H and H2 
significantly affects the observed spectrum below 1100 A. 
A direct visualization of the influence of the interstellar 
absorption on the spectrum can be seen by comparing 
similar stars in the Galaxy (Pellerin et al. 2002) and the 
Magellanic Clouds (Walborn et al. 2002). In order to 
illustrate this effect for AG Car we also display on Fig. 10 
the model spectrum taking into account the interstellar 
absorption by H and H2, using a procedure described by 
Herald et al. (2001), and assuming log[7V(H) cm^] = 21.0 
and log[iV(H2) cm^] — 21.0. It is beyond the scope of 
this work to provide a detailed analysis of the interstellar 
spectrum towards AG Car. 

6.5. Luminosity, bolometric correction, and reddening 

The luminosity of AG Car during the epochs presented 
in this paper was constrained by comparing the observed 
flux of each epoch (Table 1) with the reddened flux pre- 
dicted by the CMFGEN model (Fig. 11). For the mini- 
mum phase of 1985-1990, we used observations from the 
ultraviolet (1200 A) to the near-infrared (L band), while 
for the 2000-2003 minimum, we had available photome- 
try in the V, J, H, K, and L bands. A sole far-ultraviolet 
flux-calibrated spectrum obtained in 2001 May was ad- 
ditionally used for this epoch. Note that no AAVSO ob- 
servations were used to derive a V-band flux due to their 
relatively low accuracy compared to LTPV and ASAS-3 
photoelectric observations. 

We obtained a bolometric luminosity for AG Car of 
= 1.5 X 10^ Lq during 1985-1990 and 2000-2001, 
which is listed in Table 5 together with the other wind 
and stellar parameters. However, a lower luminosity of 
= 1.0 X 10^ Lq and = 1.1 x 10^ Lq was derived for 
the 2002 March- July and 2003 January observations, re- 
spectively. Figure 11 shows a clear trend towards lower 
luminosity during the end of the 2000-2003 minimum. 
Excluding errors in the distance, the typical error in 
Li, solely from the spectroscopic analysis is around 10% 
(gray region in Fig. 11). 

Because of our limited spectroscopic and photomet- 
ric sample, the scope of this paper is to describe the 
long-term behavior of the stellar parameters of AG Car. 
Therefore, we did not attempt to obtain a different 
bolometric luminosity for epochs relatively close in time 
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Fig. 9a. — Comparison between the observed spectrum of AG Car from 2001 April (black line) and the best model obtained with 
CMFGEN (red) in the spectral region between 3680-11000 A. For all plots, the x-axis is the wavelength in Angstroms, while the y-axis 
represents the continuum-normalized flux. The identification of the stronger spectral lines is also shown [see the electronic version of the 
journal for fits to 2002 and 2003 data]. 



(Af ~ 3 — 6 months), when httle spectroscopic variabil- 
ity was seen. This was the case, for instance, for 1989 
March-1989 December, 2001 April-2001 June, and 2002 
March-2002 July. Although AG Car clearly presented 
photometric variability during these time periods, we as- 
sumed a constant bolometric luminosity for each of these 
periods and derived an average bolometric luminosity 
which, within the errors, could describe the SED dur- 
ing the whole time interval. Note that, according to Fig. 
11, such an assumption does not change the main result 
of this Section, i.e., a lower bolometric luminosity as the 
star is moving towards maximum. 

A much larger amount of spectroscopic and photomet- 
ric data with a higher time sampling, together with a 
careful treatment of time-dependent effects in the radia- 
tive transfer model, are needed to study the behavior of 
the stellar parameters of LBVs on timescales of the or- 
der of a couple of months. Nevertheless, if the bolomet- 
ric luminosity is variable on these short timescales, its 
amplitude is likely much smaller (A log[L*/LQ] ~ 0.02 
dex) than the amplitude of the long-term variability 
(A log[L^/L0] ~ 0.17 dex). Therefore, small fluctuations 



of ±0.02 dex in \og{L-^ / Lq) may occur within each range 
of epochs shown Table 1, but such variations are within 
the errors of our modeling. 

Since we derived based on the comparison between 
the observed and model spectral energy distributions, 
different values of could in principle be obtained by 
changing M and T^s- However, unrealistic changes in 
these parameters are required in order to reach the same 
value of Li, for all the different minimum epochs; for in- 
stance. Toff would need to be reduced/increased by more 
than 2000 K for a given epoch, which is not supported 
by the line fits, in order to change the derived value of 
by 30%. A similar drastic change in M by more than 
a factor of 2, which again is not supported by the line 
fits, would be required in order to change the continuum 
emission and thus fit the data with the same for all 
epochs. Even then, very peculiar reddening laws, which 
would be different for each epoch, would be required in 
order to reproduce the shape of the SED. 

Therefore, the decrease in the bolometric luminosity 
obtained during the end of the minimum is robust, im- 
plying that the bolometric luminosity changes during the 
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Fig. 9b. — Similar to Fig. 9a, but for data taken in 2002 March. 



S-Dor cycles of AG Car. The similar maximum value of 
Li, derived for AG Car, for both minimum phases ana- 
lyzed in this work, is indicative that this is the maximum 
luminosity of AG Car during minimum. On the other 
hand, the change in bolometric luminosity is contrary 
to the currently widely accepted paradigm of S-Dor vari- 
ability at constant L^,. A detailed analysis of the changes 
of during the full S-Dor cycle will be presented in Pa- 
per III. 

The bolometric correction (BC) in the V-band for each 
epoch was also obtained from our modeling and is pre- 
sented in Table 5. As a consequence of the different stel- 
lar parameters, the BC of AG Car is also different in con- 
secutive minima, ranging from ~ —2.50 mag (1985-1990) 
to -2.15 mag (2000 July-2001 June). Further change 
in BC occurs when AG Car is moving towards maxi- 
mum {BC ~ —1.68 mag during 2002 March- July and 
BC ~ -1.22 mag in 2003 January). 

While it is well known that the bolometric correction 
of stars which have a negligible wind (such as the Sun 
and OB stars) is dependent on Tcs and logg, we found 
that the bolometric correction of LBVs, because of the 
atmospheric extension, is dependent not only on Tcs and 
logg but also on the wind density (and thus on M and/or 



Woo)- Since Tcff also depends on the wind density, LBVs 
with different stellar parameters can have a similar value 
of Tcff but significantly different BCs; this is exactly the 
case for AG Car during 2001 April and 2002 March (see 
Table 5). While T^a was close to - 17,000 K during 
both epochs, the star was very different at depth, as can 
be readily evidenced by the very different temperatures 
at the hydrostatic radius in 2001 April (T* ~ 21, 900 K) 
and 2002 March (T* ~ 18,700 K). Therefore, the BCs 
were different, and the effective temperature was similar 
just by coincidence. A complete analysis of the bolo- 
metric correction of LBVs over a broad range of physical 
parameters is beyond the scope of this work, and will be 
presented elsewhere. 

Together with L*, the reddening law towards AG Car 
was obtained through the comparison of the observed 
flux and CMFGEN models reddened according to the 
law of Fitzpatrick (1999). The best fits were obtained 
using E{B -V) = 0.65 ± 0.01 and Ry = 3.5 ± 0.1, and 
are shown in Figure 11. Taking the uncertainties into ac- 
count, it was possible to fit the data from both minima 
with the same reddening law, corroborating the physi- 
cal parameters obtained for both epochs. The value of 
E{B — y) = 0.65 is only slightly higher than previous 
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Fig. 9c. — Similar to Fig. 9a, but for data taken in 2003 January. Note that Ha is saturated and was not used as a diagnostic. 



determinations {E{B — V) = 0.63) obtained using sim- 
ple models and assuming constant bolometric luminosity 
(Humphreys et al. 1989; Shore et al. 1996). 

Although the value of Ry = 3.5 obtained in this work 
is higher than the standard value due to diffuse dust 
= 3.1), it is lower than Ry = 3.9, as suggested by 
Schulte-Ladbeck et al. (1994) and based on the analy- 
sis of the polarimetric data of AG Car. Our results are 
compatible with the presence of large grains (~ 1/im) 
in the circumstellar environment of AG Car (Hyland & 
Robinson 1991), and with the very low amount of nebular 
reddening which has been inferred {Ay ^ 0.025 — 0.05, 
McGregor et al. 1988a). 

The agreement between the CMFGEN models and the 
observed SED of AG Car from the far-ultraviolet to the 
L-band is superb. The free-free emission of the wind 
predicted by the CMFGEN models agrees very well with 
the observed SED in the near-IR, suggesting that the 
amount of warm and hot dust 800-1500 K) present 
in the close vicinity of AG Car (1-2 ") is negligible. 

7. DISCUSSION 

7.1. Chemical abundances on the surface of AG Car 
and on the bipolar nebula 



The chemical abundance on the surface of AG Car, 
presented in Table 4, clearly shows enhancement of He, 
N, and Na, together with a depletion of H, C, and O 
compared to the solar value, which is a fingerprint that 
CNO-processed material is present on the stellar surface 
and wind of AG Car. Since AG Car is a prototype of the 
LBV class, the abundances obtained in this work quanti- 
tatively confirm the prediction that LBVs have material 
processed by the CNO-cycle and significant He enrich- 
ment on their surfaces (see e.g. Humphreys & Davidson 
1994, Crowther 1997), which put tight constraints on 
their evolution. 

Table 6 presents the abundances obtained in this work 
for the surface of AG Car and the available nebular val- 
ues. Their comparison might provide clues to the evolu- 
tion of AG Car during the LBV phase, since the nebula 
was ejected ^ lO"' years ago (Smith et al. 1997; Earners 
et al. 2001). Unfortunately, nebular abundances are only 
reliable for N and O, since the He and C lines are weak or 
absent in the nebular spectrum (Mitra & Dufour 1990; 
de Freitas Pacheco et al. 1992; Nota et al. 1992; Smith 
et al. 1997). The He abundance was estimated by Earners 
et al. (2001) based on several evolutionary assumptions 
and on the ratio N/0 obtained for the nebula. 



On the nature of AG Car I. Minimum phases 



17 



0.6 
0.4 
0.2 

0.0 
0.6 

0.4 

0.2 

0.0 



AG Car FUSE obs. 2001 May 

CMFGEN model no IS abs. 

CMFGEN model + lS abs.(N„=10" cm-',N„=10" cm-') 



= 10" cm"*) ;, ;■, 





AG Car FUSE obs. 2001 May 
Fe II CMFGEN model 
Fe III CMFGEN model 



940 



960 



980 



1000 



1020 




1020 



1040 



1060 



1080 



1100 






1100 



1120 



1140 
Wavelength [A ] 



1160 



1180 



Fig. 10. — Far-ultraviolet spectrum of AG Car observed with FUSE in 2001 May (black solid hue) compared with the best CMFGEN 
model for this epoch (green dotted line), reddened using E{B — V) = 0.65 and Ry = 3.5. The CMFGEN model was scaled in flux in 
order to better match the continuum level. The best model, taking into account an interstellar absorption spectrum due to H I with 
log[Af(H) cm^] = 21.0 and H2 with log[Af(H2) cm ] = 21.0, is shown in order to infer the influence of those transitions on the morphology of 
the far-ultraviolet spectrum of AG Car (purple solid line). For each wavelength range, the bottom panel shows the best CMFGEN model 
spectrum computed including continuum emission and bound-bound transitions due to only Fe II (red solid line) or only Fe III (blue solid 
line), in order to facilitate line identiflcation. These ions are the main contributors to the line blanketing in the far-UV and illustrate how 
well CMFGEN can handle the effects caused by the extreme blending due to tens of spectral lines. 



TABLE 5 

Stellar and wind parameters of AG Car during minimum epochs 



Epoch 


ya 


BC 


log(L,/Lo)'= 




.Rphot 


T» 


T„ff 


M 






1 






(mag) 


(mag) 


(Re) 


(flo) 


(K) 


(K) 


(Moyr- 


(kms^^) 




1985 Jul - 1986 Jun 


7.98 


-2.50 


6.17 


58.5 


78.7 


26,450 


22,800 


1.9 X 10 


-5 


300 


0.10 


1987 Jan - 1990 Jun 


8.00 


-2.52 


6.17 


59.6 


78.7 


26,200 


22,800 


1.5 X 10 


-5 


300 


0.10 


1990 Dec - 1991 Jan 


7.71 


-2.23 


6.17 


67.4 


88.5 


24,640 


21,500 


1.5 X 10 


-5 


300 


0.10 


2000 Jul-2001 Jun 


7.63 


-2.15 


6.17 


85.3 


141.6 


21,900 


17,000 


3.7 X 10 


-5 


105 


0.15 


2002 Mar-2002 Jul 


7.60 


-1.68 


6.00 


95.5 


124.2 


18,700 


16,400 


4.7 X 10 


-5 


195 


0.25 


2002 Nov 


7.20 


-1.28 


6.00 


120.4 


170.4 


16,650 


14,000 


6.0 X 10 


-5 


200 


0.25 


2003 Jan 


7.03 


-1.22 


6.04 


115.2 


171.3 


17,420 


14,300 


6.0 X 10 


-5 


150 


0.25 



^Notc that the V magnitude of AG Gar is variable within each range 
of dates (see Table 1 for values for a specific date). 

^The photometric variability inside each range of epochs, if due to 
changes in L^, , imply variations in the value of log (L* / Lq ) listed above 
inside each range of epochs by ±0.02 dex. 

We support the conclusion of Lamers et al. (2001) that 
the He abundance on the surface of AG Car is higher 
than in the nebula^^. In addition, we found, despite the 
uncertainties, that the ratio N/0 is significantly higher 
on the stellar surface (N/0=39l^g by number) than in 
the nebula (N/0=6± 2, Smith et al. 1997). Therefore, 
we conclude that the surface of the star has more CNO- 
processed material than the nebula. This is consistent 
with the idea that the massive nebula is composed of 
less chemically-processed outer layers which were ejected 
from the central star. The pre-existing interstellar bub- 

Lamers et al. (2001) used the He abundance of AG Car derived 
by Smith et al. (1994), which is very similar to the value obtained 
in this work. 



ble likely contains negUgible amount of material com- 
pared to the total mass of the nebula (Lamers et al. 
2001). Whether all the nebular mass was ejected in a 
single giant eruption ~ 10* years ago relics on the knowl- 
edge of the mass-loss history of AG Car earlier than and 
during the LBV phase, which is unknown. 

The N/0 ratio obtained for AG Car is roughly in 
line with the predictions of the evolutionary models by 
Meynet & Maeder (2003). A star with an initial mass 
of 85 Mq reaches a similar He content as AG Car at an 
age of about 3x10^ years, when the evolutionary models 
predict N/0~ 26 and N/0~ 20 for the non-rotating and 
rotating cases, respectively. Detailed comparisons be- 
tween AG Car and evolutionary tracks, and the implica- 
tions of our results for the current evolutionary scenario 
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Fig. 11. — Comparison between the observed spectral energy distribution of AG Car obtained from photoelectric data during minimum 
epochs (red asterisks) and CMFGEN models with different bolometric luminosities. For 2000 July we display the visual magnitude (red 
open diamond). The panels are ordered chronologically, and in decreasing effective temperature, from left to right, top to bottom. For each 
panel we show selected models with = 2.2 X 10® Lq (dark green dash-triple dotted line), L* = 1.5 X 10® Lq (black dash-dotted line), 
Li, = 1.1 X 10® Lq (Ught green sohd hne), = 1.0 X 10® Lq (blue solid line), and = 0.7 X 10® Lq (purple dashed hne). For a given 
epoch, models with different had their ij* and M scaled in order to still fit the observed spectral lines. The gray region corresponds to 
the uncertainty in the model flux which comes from the luminosity. In all panels, the CMFGEN models were reddened using the Galactic 
extinction law of Fitzpatrick (1999) with Ry = 3.5 and E{B — V) = 0.65. For each epoch, there is superb agreement between the best 
model and the observed data using the same reddening law, although different luminosities are needed. Note the trend towards lower 
luminosity as the effective temperature (radius) decreases (increases) during 2001-2003. The uncertainty in the observed photometry is 
smaller than the size of the symbols in the UV, optical, and in the near-IR. 



of massive stars will be presented in Paper II. 

In the last lO'^ years, additional CNO-processed mate- 
rial could have been brought to the surface of AG Car. 
The star is likely a fast rotator (Groh et al. 2006), sug- 
gesting that rotational mixing could efficiently occur due 
to meridional circulation and shear instability (Maeder 
& Meynet 2000). In addition, approximately 0.5 Mq 
has been removed from the star due to the steady stellar 
wind since the nebula was ejected, assuming an average 
mass- loss rate of 5 x 10~^ M0yr~^ during the last 10'' 
years (Paper III). 

7.2. Comparison with previous works 

Fundamental parameters of AG Car during the mini- 
mum of 2000-2003 have been presented here for the first 
time and, therefore, only our results for the minimum 



TABLE 6 

Comparison between surface and nebular abundances of AG Car 





He/H 


N/O 


Reference 




(number) 


(number) 




Surface 


0.43 ±0.08 




this work 


Nebula 


~0.3^ 


6±2 


Smith ct al. (1997) 



^Indirectly obtained by Lamers et al. (2001). 



of 1985-1990 can be compared with previous works. We 
summarize in Table 7 the results obtained by different 
works. 

In general, we derived a higher L*, higher T^,, lower 
Teff, and lower M compared to previous works. These 
different results are directly or indirectly related to the 
inclusion of wind clumping, full line blanketing, detailed 
treatment in non-LTE of the populations and of the ra- 
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diative field in the co-moving frame, and improved treat- 
ment of the subsonic part of the wind in our analysis. 
The lower mass-loss rates derived by us are a direct con- 
sequence of the relatively high degree of clumping found 
in the wind of AG Car (/ ~ 0.1), and the assumption 
of different L and Voo by previous authors also affects 
the value of M . The lower M and higher L^, derived by 
our work implies a reduction in the wind performance 
number rj = Mvacc/L. During the minimum of 1985- 
1990 we found a quite modest value of 77 ~ 0.16. The 
increase in in our results can be explained by the 
back-warming of the base of the wind due to the inclu- 
sion of line blanketing in our modeling. The significant 
differences in T^s also indicate that obtaining T^a from 
simple empirical relationships, such as using the V-band 
flux, do not provide a correct answer for hot luminous 
stars with dense winds, such as AG Car. 

The He abundance obtained in this work (He/H= 
0.43 ± 0.08) is compatible with the value of He/H ^ 0.41 
determined by Smith et al. (1994). All the other abun- 
dances have been derived here for the first time, including 
the detection of CNO-processed material on the surface 
of AG Car. 

The wind terminal velocity derived by us (woo — 
300 kms~^) is higher than the value of v^a — 250 kms~^ 
obtained by previous works using the FWHM of optical 
fines (Leitherer et al. 1994; Smith et al. 1994; Stahl ct al. 
2001), because wc used the saturated absorption profile 
of resonance lines in the UV to infer v^o , which is a more 
accurate diagnostic. We found that the optical lines are 
sensitive to the velocity law, and usually provide a lower 
limit to Woo- Guo & Li (2007) derived (3 = 0.5 for AG Car 
during 1990 based on the fitting of the continuum flux, 
but under the assumption of a series of fixed parameters 
(M, i?*, Teff, L*, Voo) obtained from Stahl et al. (2001). 
All these parameters, as well as wind clumping and line 
blanketing, affect the shape of the continuum, and the 
assumed values by Guo & Li (2007) are significantly dif- 
ferent than the values derived in this paper. Our detailed 
fits to the line profiles require a much slower velocity law 
(/3 = 3) for the wind of AG Car than that determined by 
Guo & Li (2007). 

7.3. On the nature of the changes in the bolometric 
luminosity during the S-Dor cycle 

The evolution of LBVs during the S-Dor cycles is 

widely assumed to occur roughly at constant bolomet- 
ric luminosity (van Gendercn 1979; Stahl et al. 1983; 
Humphreys et al. 1989; Leitherer et al. 1989, 1994; Smith 
et al. 1994; Shore et al. 1996; de Koter et al. 1996; 
Stahl et al. 2001; Vink & de Koter 2002; Walborn et al. 
2008). However, the few previous works that investigated 
changes in the bolometric luminosity of LBVs during the 
S-Dor cycle (e.g. Humphreys et al. 1989 for AG Car, 
and Walborn et al. 2008 for R127) employed standard 
bolometric corrections of OB stars, which we argue are 
not realistic for LBVs, and photometry restricted to the 
F-band. For instance, the bolometric corrections ob- 
tained in this work for AG Car arc more negative by 
0.15 to 0.60 mag, depending on the epoch, than those 
of normal B supergiants (with a similar low logg) ob- 
tained using line-blanketed models such as CMFGEN 
(Crowther et al. 2006) or TLUSTY (Lanz & Hubeny 



2007). When compared to the BCs of B supergiants pro- 
vided by works which employed non-blanketed models, 
such as Humphreys & McElroy (1984), the BCs derived 
for AG Car are even more discrepant, being more nega- 
tive by 0.40 to 0.90 mag, depending on the epoch. 

Through a detailed multi-wavelength spectroscopic 
modeling using CMFGEN, we foimd in Sect. 6.5 that the 
bolometric luminosity of AG Car decreases by a factor 
of 1.5 (~ 0.2 dex) from minimum towards the maximum 
phase of the S-Dor cycle Sect. 6.5. Such a finding de- 
serves further discussion, in particular because the vari- 
able Li, could provide a crucial constraint on the mecha- 
nism driving the S-Dor type variability. 

A reduction of the bolometric luminosity during the 
S-Dor cycle was also obtained by Lamers (1995) for the 
LBV S Doradus, and the result was interpreted in terms 
of the energy used to expand the outer layers of the star 
from minimum to maximum. Lamers (1995) found that 
about 10^'^ to 10^^ Mi, are involved in the expansion 
of S Doradus, corresponding to ~ 0.17 Mq (assuming 
~ 45 Mq, Lamers 1995). 

We propose that a similar behavior is present in 
AG Car. In order to obtain an order-of-magnitude es- 
timate of the amount of mass taking part in the ex- 
pansion of the star, we assume that L^, decreased from 
Li, = 1.5 X lO'^ Lq to Li, = 1.0 X 10^ Lq on 2001 June 15 
(the last epoch for which we derived = 1.5 x 10^ Lq, 
see Sect. 6.5), and remained constant until 2003 January 
11. The elapsed time At between 2001 June 15 and 2003 
January 11 is 575 days and, therefore, the missing radia- 
tive energy during these epochs is Ai?rad = AL ■ At ~ 
1 X lO"'^ erg. Assuming that the layer is expanding from 
85 Rq to 115 Rq, and following Lamers (1995), 

^^'^"-^R^ \85-TT5/' 

where Mcs = M*(l - F). Since our CMFGEN models 
predict that F ~ 0.8, we obtain M^g ~ 14 Mq assum- 
ing Mi, ~ 70 Mq (Paper II). Using Eq. 2, we derive 
that roughly 8.3 x 10"^ ~ 0.6 Mq is involved in the 
expansion of the star during the S Dor cycle. 

This result should be viewed with caution, since an 
uncertainty of at least a factor of 3 is present due to 
the various assumptions made above. The amount of 
~ 0.6 Mq can be easily increased by a factor of 2 if we 
include the thermal energy released by the expanding 
layer (sec discussion in Lamers 1995). If we also take 
into account that wc did not analyze in this paper the 
epochs corresponding to the maximum in the lightcurve 
of AG Car, which was only reached in 2004-2005 (Paper 
III), an additional amount of energy would be required, 
increasing by another factor of ^ 2 the amount of mass 
involved in the expansion. Therefore, although detailed 
modeling is obviously needed, it might be possible that 
0.6 to 2 Mq are taking part in the expansion of AG Car 
during the S-Dor cycle. 

Interestingly, such an amount of mass is an order of 
magnitude lower than the nebular mass found around 
AG Car (- 15 - 30 Mq, Voors et al. 2000) and that of 
the Homunculus nebula around Eta Car (^ 12 — 20 Mq, 
Smith et al. 2003). On the other hand, the amount of 
mass involved in the S-Dor type variability of AG Car 
(this work) and S Doradus (Lamers 1995) is comparable 
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TABLE 7 

Fundamental parameters of AG Car obtained by different works 



1989 March 

Work log L*/Lq T\ M/-"-5 ^ 7 Hc/H Vj 

(K) (K) (10~^ M0yr~l) (kms~l) (number) (mag) 

this 637 26,200 22,800 4?7 300 Ol 043 OIT -8.17 
SCP94^ 6.04 26,000 5^6 250 1.0 0.41 0.60 -7.7 



1990 December 

Work log L^/Lq n r,ft M/-0-5 / He/H v My 

(K) (K) (10~^ Mgyr-i) (kms'^) (number) (mag) 

this 6.17 24,640 21,500 4.7 300 0.1 0.43 0.17 -8.46 
LEI94'' 6.00 ... 21,000 10.0 250 1.0 1.00 1.42 -8.0 

STAOl'^ 6.00 24,060 3^3 225 1.0 0.41 0.35 -8.0 



^Smith ct al. (1994). 
''Leithercr ot al. (1994). 
<=Stahl ot al. (2001). 

to the nebular mass found around many low-luminosity 
LBVs and to that of the Little Homuneulus around Eta 
Car (see Smith & Owocki 2006 and references therein 
for a compilation of the amount of nebular mass found 
around LBVs). Such similarity might suggest a link be- 
tween some of the giant outbursts of LBVs and their 
S-Dor cycles, and wc speculate that the S-Dor cycles 
could actually be failed Giant Eruptions. Instead of los- 
ing several solar masses like Eta Car did in the 1840's, 
we suggest that during the S-Dor type variability such 
an amount of mass is never ejected from the star. 

8. CONCLUSIONS 

The detailed spectroscopic analysis of AG Car using 
the radiative transfer code CMFGEN has provided ad- 
ditional insights on the LBV phenomenon and on the 
S-Dor type variability. Below we summarize the main 
conclusions of this paper. 

1. Following the determination of a high content of 

He, N, and Na, and the depletion of H, C, and 
O, compared to the solar values, we infer that 
CNO-proccsscd material is present on the surface of 
AG Car. The actual surface is probably the peeled- 
off inner layers of the star, which are exposed due 
to the action of a continuous stellar wind and the 
ejection of the outer layers due to a giant erup- 
tion that occurred 10** years ago. Since AG Car is 
likely a fast rotator (Groh et al. 2006), rotational 
mixing might also have played a role to increase the 
content of CNO-processed material on the surface. 
The abundances obtained for Si, Al, Mg, Fe, Ni, 
Co, Cr, and Mn are consistent with a solar abun- 
dance. 

2. The minimum phases of the S Dor cycle of AG Car 
are not equal to each other. The consecutive mini- 
mum phases of 1985-1990 and 2000-2003 were dif- 
ferent in duration, visual magnitude, stellar tem- 
perature, mass-loss rate, and wind terminal veloc- 
ity. We suggest that these differences arise due 
to different underlying stellar parameters, which 
cause different mass-loss rates and wind terminal 
velocities. In Paper II we will analyze whether 
the bi-stability mechanism (Pauldrach & Puis 1990; 
Vink & de Koter 2002) can explain this behavior. 

3. The duration of the last two minimum phases of the 



5 Dor cycle of AG Car and their variability are re- 
lated to the maximum temperature achieved. The 
hotter minimum of 1985-1990 was roughly twice as 
long as the cooler minimum of 2000-2003. On the 

other hand, AG Car was markedly more variable 
during 2000 2003 than during 1985 1990. 

4. The comparison between the observed flux of 
AG Car from the ultraviolet to the near-infrared 
with the model flux yielded a color excess of E{B — 
V) = 0.65±0.01, with Ry = 3.5±0.1. The value of 
Rv is consistent with the presence of large grains in 
the circumstellar environment of AG Car (Hyland 

6 Robinson 1991), but it is lower than the value 
of Rv = 3.9 proposed by Schulte-Ladbeck et al. 
(1994). 

5. The maximum value of the bolometric luminos- 
ity obtained during minimum phases was L^, = 
1.5 X IO^Lq for an assumed distance of d = 6 kpc. 
Excluding the uncertainty in the distance, the er- 
ror in due to the spectroscopic analysis is around 
10%. 

6. Contrary to the current paradigm, we found that 
Li, decreases by a factor of 1.5 from minimum as 
the star moves towards visual maximum of the S- 
Dor cycle. Since further reduction in Teff occurs 
during maximum (Paper III), it might be possible 
that decreases even more during the maximum 
of the lightcurve. 

7. Assuming that the decrease in the bolometric lu- 
minosity of AG Car is due the energy used to ex- 
pand the outer layers of the star (Lamers 1995), 
we found that roTighly 0.6 — 2 Mq are involved in 
the expansion. Although that is much lower than 
the nebular mass found around AG Car and that 
the mass of the Homuneulus around Eta Car, the 
amount of mass involved in the S-Dor type vari- 
ability of AG Car is comparable to the nebular 
mass found around many low-luminosity LBVs and 
to that of the Little Homuneulus around Eta Car. 
Such similarity might suggest a link between some 
of the giant outbursts of LBVs and their strong 
variability during the S-Dor cycles. We speculate 
that the S-Dor type instability could be a failed Gi- 
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ant Eruption, with the several solar masses never 
being released from the star. 
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